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In the present work we propose to study neutrino oscillations employing sources of monoenergetic 
neutrinos following electron capture by the nucleus. Since the neutrino energy is very low the smaller 
of the two oscillation lengths, L23, appearing in this electronic neutrino disappearance experiment 
can be so small that the full oscillation can take place inside the detector and one may determine 
very accurately the neutrino oscillation parameters. Since in this case the oscillation probability 
is proportional to sin 2 26n, one can measure or set a better limit on the unknown parameter 613. 
This is quite important, since, if this mixing angle vanishes, there is not going to be CP violation in 
the leptonic sector. The best way to detect it is by measuring electron recoils in neutrino-electron 
scattering. One, however, has to pay the price that the expected counting rates are very small. 
Thus one needs a very intensive neutrino source and a large detector with as low as possible energy 
threshold and high energy and position resolution. Both spherical gaseous and cylindrical liquid 
detectors are studied. Different source candidates are considered. 

PACS numbers: 13.15. +g, 14.60Lm, 14.60Bq, 23.40.-s, 29.25.Rm, 95.55.Vj. 

INTRODUCTION. 

The discovery of neutrino oscillations can be considered as one of the greatest triumphs of modern 
physics. It began with atmospheric neutrino oscillations [l[ interpreted as — > v r oscillations, as 
well as v e disappearance in solar neutrinos These results have been recently confirmed by the 
KamLAND experiment 0], which exhibits evidence for reactor antineutrino disappearance. As a 
result of these experiments we have a pretty good idea of the neutrino mixing matrix and the two 
independent quantities Am 2 , e.g |m| — m 2 | and |m§ — m||. Fortunately these two Am 2 values are 
vastly different, 

Am 2 ! = |m 2 - m\\ = (7.65l°; 23 ) x 

and 

Am 2 2 = |m 2 - m 2 | = (2.4±°;i 2 ) x W\eV) 2 . 

This means that the relevant L/E parameters are very different. Thus for a given energy the 
experimental results can approximately be described as two generation oscillations. For an accurate 
description of the data, however, a three generation analysis (3] , 0] is necessary. 

In all of these analyses the oscillation length is much larger than the size of the detector. So one 
is able to see the effect, if the detector is placed in the right distance from the source. It is, however, 
possible to design an experiment with an oscillation length of the order of the size of the detector. 
In this case one can start with zero oscillation near the source, proceed to maximum oscillation near 
the middle of the detector and end up again with no oscillation on the other end. This is achieved, 
if one considers a neutrino source with as low as practical neutrino energy. The main requirements 
are as follows: 

• The neutrinos should have as low as possible energy so that the oscillation length can be 
minimized. At the same time it should not be too low, so that the neutrino-electron cross 
section is sizable. 

• A monoenegetic source has the advantage that some of the features of the oscillation patterns 
are not washed out by the averaging over a continuous neutrino spectrum. 
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• The life time of the source should be suitable for the experiment to be performed. If it is 
too short, the time available will not be adequate for the execution of the experiment. If it is 
too long, the number of counts during the data taking will be too small. Then one will face 
formidable backgrounds and/or large experimental uncertainties. 

• The source should be cheaply available in large quantities. 

Clearly a compromise has to be made in the selection of the source. It appears, unfortunately, that 
the oscillation length one may have to live with, may be quite a bit longer than that ofthe previously 
considered triton source with a maximum energy of 18.6 keV, i.e. an average oscillation length is 
6.5m. 

At such low energies the only neutrino detector, which is sensitive to neutrino oscillations, is one, 
which is capable of detecting recoiling electrons. With this detector one is sensitive to all neutrino 
flavors I n fact one can detect: 

• electrons which are produced by electron neutrinos via both the charged and the neutral current 
interaction. These will manifest as electron neutrino disappearance 

• electrons are produced from the other two neutrino flavors due to the neutral current interac- 
tion. These two flavors are due to the muon and tau neutrino appearance. 

Since the importance of the cross section due to the charged current relative to that due to the 
neutral current depends on the electron energy, one has a novel feature, i.e. the effective oscillation 
probability depends on the electron energy. Thus the results may appear as disappearance oscillation 
in some kinematical regime and as appearance oscillation in some other regime. In this paper we will 
examine how these features can best be exploited in determining the neutrino oscillation parameters. 



ELASTIC NEUTRINO ELECTRON SCATTERING 

For low energy neutrinos the historic process neutrino-electron scattering 0] is very useful. 
The differential cross section [8[ takes the form 

da ( da\ ( da \ 

The second term, due to the neutrino magnetic moment, is inversely proportional to the electron 
energy and, at the low electron energies of the present set up, may be used in improving the current 
limit of the neutrino magnetic moment by two orders of magnitude. It is not, however, the subject 
of the present study, which is concerned with neutrino oscillations. The cross section in the rest 
frame of the initial electron due to weak interaction alone becomes: 



%) =^"<> "" 2 



(9v--w) 2 [l-^] 2 + (9j-9f-) ! |fl 



gv = 2 sin 2 6w + 1/2 (v e ) , gy = 2 sin 2 9 W — 1/2 {v^,v T ) 



g A = \/2 (v e ) , g A = -l/2 (va,v T ) 



For antineutrinos g A — > ~9a- 

The scale is set by the week interaction: 



'-F^ = 4.45 x lfj- 48 — (3) 
2tt keV y ' 
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The above equation must be modified, if one wishes to include the fact that the initial electron is 
bound. For neutrino energies in the tens of keV, however, it has been shown by Gounaris, Paschos 
and Porfyriadis @ that such a modification lowers the cross section by no more than 10%. So we 
are not going to consider such effects here. 

The electron energy depends on the neutrino energy and the scattering angle and is given by: 



T = 



2 m e {E v cos Of 



{m e + E v f - (£ !/ cos0) 2 
The last equation for sufficiently low energies can be simplified as follows: 

2(£ i/ cos0) 2 



(4) 



(5) 



The maximum electron energy depends on the neutrino energy squared. The total cross section as 
a function of neutrino energy is shown in Fig. [T] 
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FIG. 1: The total (v e ,e ) cross section in the absence of oscillations as a function of the neutrino energy. 
In the present work the neutrinos have a definite energy. 



NEUTRINO SHORT BASELINE OSCILLATIONS 

The electron neutrino, produced in weak interactions, can be expressed in terms of the three mass 
eigenstates as follows: 

v e = cos 6*i2 cos 6 13 v\ + sin 6*12 cos 6*13 V2 + sin #13 e lS V3 (6) 
Then the v e disappearance oscillation probability is given by: 

P(u e -> v e ) = 1 - [cos 4 0i3 sin 2 26»i2 sin 2 (tt-^-) 

+ sin 2 0i2 sin 2 20 i3 sin 2 (tt-—) 

L32 

+ cos 2 0i2 sin 2 20i3 sin 2 (tt— )] (7) 

^31 
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with 

^ij = 2 2 (^) 

mf — ml 

sinOis is a small quantity constrained by the CHOOZ experiment and it will be further investigated 
below. In the limit in which it is neglected, one recovers the standard expression familiar from the 
solar neutrino oscillation analysis, if we take s\2 ~ sin9 so i ar , c\i ~ cos9 so i ar , S oiar as determined 

31 ~ ^ m 32 



from the solar neutrino data. Assuming Amjj m Am| 2 we find 



P{u e -> v e ) rj 1- [(sin20 soiar ) 2 sin 2 (7r-^) + 

-^21 

sin 2 (20 13 )sin 2 (7r— )] (9) 

^32 

where the term proportional to sin 2 (2^13), connected with the small oscillation length [10] L32, 
(L 2 i ~ 32L 32 ), is relevant for the experimental approaches discussed in the present paper. As we 
have already mentioned in connection with the NOSTOS experiment [lO], in addition to the electronic 
neutrinos, which are depleted by the oscillations, is sensitive to the other two neutrino flavors, which 
can also produce electrons via the neutral current interaction. These flavors are generated via the 
appearance oscillation: 

P{ve -> ^ v a ) w (sin26> so ; ar ) 2 sin 2 (7r-^) + 

sin 2 (20 13 )sin 2 (^-^) (10) 

-^32 

Thus the number of the scattered electrons, which bear this rather unusual oscillation pattern, is 
proportional to the (v e ,e~) scattering cross section with a proportionality constant C{E Vl T) given 
by: 

C{E V ,T) = l- X {E v ,T) 

(sin26U Zar ) 2 sin 2 (7r-^) + sin 2 (26> 13 ) sin 2 (7r-^). 

-^21 -^32 

(11) 

The function x(.E„,T), which represents the relative difference between the cross section of the 
electronic neutrino and its other two flavors, has been previously discussed [l0| . 



GEOMETRIC CONSIDERATIONS 

The total neutrino electron scattering cross section as a function of x and L can be cast in the 
form: 

a(L,x) = a(0,x)(l- X (x) P (L,x)) (12) 

with x = — and 

m e 

, , Glm 2 „x 2 (l7.7464a: 2 + 15. 3098a; + 3.36245) , , 

^ X) = ^r— (2XTW (13) 



is the total cross section in the absence of oscillations. Furthermore 

I595922L\ 
330a; J 
595922L\ 
"lCte J 



, 1 0.0595922L\ . 2 . n 
p(L,x) = sm [ Ism (20 so ; ar ) + 

, l 0.0595922L\ . , ,„„ . 
sm 2 sm 2 (20i 3 ) (14) 
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with L the source detector distance in meters and 

2.8664a; 2 + 4.1498a; + 1.50245 
X ^ X ' ~ 17.7464a; 2 + 15.3098a; + 3.36245 ( ' 

We will assume that the volume of the source is much smaller than the volume of the detector. 
We will consider two possibilities: 



Spherical detector with the source at the origin 



The number of events between L and L + dL is given by: 



dN = N v n e 4W A L ' dL a(L, x) = N v n e dLa(L, x) (16) 



or 

dN 



= N v n e (j(L 1 x) (17) 



To compare with other geometries we rewrite this as follows: 

dN 



. = R N u n e cr(x, p) (18) 
dp 



or 



= -^R N„n e9s (p)a{x,p), g s ( P ) = l (19) 

with N v the number of neutrinos emitted by the source, n e the density of electrons in the target, 
Rq the radius of the target and p = L/R . a(x, p) is the neutrino - electron cross section in units of 
G 2 F m1/2ir. The geometry factor g s (p) is in this case unity. 

Cylindrical detector with the source at the center of one of its basis 

The number of events between r and r + dr and z and z + dz is now given by: 

dN = N v n e 2 . a(x,Vr 2 + z 2 ) 20 

ATT(r z + z z ) 



which yields 



~f~77 = N v n e R\ l\ 2 a(x, -VC 2 +«V, (21) 
dpdQ 2 Q z + u z p z u 



where R is the radius of the cylinder, u = R/h (h is the length of the cylinder), p = r/R and 
( = z/h. This can be written as 

dN , , R 



= Ag c (p, C, u, R)-a(x, - ^C 2 + u z p z ), (22) 



dpdC, ' 2 u 



g c (p,(,u,R)= 2 , K=±J^RN v n e (23) 

Q z + u z p z 2ir 

Note that the geometric factor g c {p, (, u, R), which is absent in the spherical geometry, in this case 
is less than unity. 
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It is instructive to compare the two geometries in the absence of oscillations. In such a case we 
find that the number of events is given by: 



, J / 1 ( 1 A log (u 2 + 1) \ 
iV c = N v n e R<r{x^)- I 2tan~ 1 f - J + V '- \ 



(24) 



while for a half sphere of the same volume we find: 



N s = N u n e Ra(x, 0) 



3/3 

V 2 

2?/u 



(25) 



The relative merit N c /N s is shown in Fig. [5] This ratio becomes unity for u — 1.12, which is much 
larger than the value of 11/91 of the planned detector. Note, however, that for u=ll/91 the relative 
merit is 0.75, .i.e. not far from unity. 
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FIG. 2: The relative geometric merit, N c /N s 
volume as a function of u = R/h 



of a cylinder compared to that of a half sphere of the same 



CANDIDATES FOR NEUTRINO OSCILLOMETRY 



One can ask whether the relevant candidates for small length oscillation measurements exist in 
reality. The analysis of known information on nuclides show that there are many cases with a small 
neutrino energy in the orbital electron capture by nucleus. Since this process has the two-body 
mechanism, the total neutrino energy is equal to the difference of the total capture energy QeC and 
binding energy of captured electron Bi : 

Qe = Qec ~ Bi (26) 
This value can be easily determined because the capture energies are usually known (or can be 



measured very precisely by the ion-trap spectrometry ll[) and the electron binding energies are 
tabulated |12j . The main feature of the electron capture process is the monochromaticity of neu- 
trino. This paves the way for the neutrino oscillometry, which is introduced below. Table U shows 
the relevant nuclides, which emit monoenrgetic neutrinos with energy less than 750 keV, half-life 
not more than 10 5 y, and whose production rates in the neutron reactors can provide rather high 
intensity of neutrinos (> 10 12 s _1 ). Columns 2-4 of Table U show the decay characteristics of the cor- 
responding nuclides, column 5 gives a half value of the small oscillation length L32, column 6 shows 



7 



the recoil electron maximal energy in the neutrino scattering on electron, and column 7 presents the 
masses of samples that can reasonably be produced by irradiation of natural targets in reactors with 
the assumed neutron flux of 5 x 10 14 cm _2 s _1 and known neutron capture cross-sections for individ- 
ual nuclides. We can see from Table Q] that there is a variety of nuclides with different half-lives and 
capture energies. Short half-life nuclides, as e.g. 51 Cr, 55 Fe, 71 Ge, and 103 Pd, have an advantage in 
the sense that they can be rapidly produced in reactors, whereas long-lived ones, as 41 Ca, 157 Tb, 
163 Ho and 193 Pt need long-term irradiation in order to accumulate necessary decay intensity. The 
nuclides 157 Tb and 163 Ho can be accumulated after disintegration of 157 Dy and 163 Er, respectively, 
which in turn can be produced by neutron irradiation of natural dysprosium and erbium or enriched 
samples of 156 Dy and 162 Er. To produce a high intensity sources one needs a huge amount of ir- 
radiated material. For the production of 71 Ge with neutrino intensity of 2 x 10 18 s _1 (55 MCi) the 
300 kg of natural Ge should be used with the irradiation time of 20 days at the reactor. Obviously, 
these irradiations in principle can be repeated many times. For the production of 51 Cr ca. 40 kg 
enriched target of 50 Cr used for GALLEX-detector calibration [l3| can be favorably borrowed. The 
long half-life implies the production of smaller source intensity even after a long-time neutron irra- 
diation, which can not be fully compensated by subsequent long-time data acquisition. Meanwhile, 
just these long-lived nuclides with the corresponding small neutrino energies can provide as small as 
possible short oscillation lengths L 32 . 

The values of L32 in the column 5 were determined by the formula of Eq. ([8]) . Thus we can write: 

L ^ = ^Smr^ L32H ^ E4keV] (27) 

The values in the square brackets in Eqs (|2"T|) indicate the dimensions used. As can be seen from 
Table U (column 5) one can roughly divide the nuclides presented there into two categories: those 
which have L32 < 50 m and those with L32 > 110 m. For the former nuclides the TPC counting 
method can be used in the gas-filled NOSTOS sphere approach (lol |. 14 1, whereas for both and 



mainly for the latter category with the larger L32, the long liquid scintillator (LS) detector [15| is 
more preferable. Both methods are discussed in a forthcoming section. The goal of both approaches 
is to scan the monoenergetic neutrino-electron scattering events by measuring the electron recoil 
counts in a function of distance from the neutrino source prepared in advance at the reactor/s. 
This scan means point-by-point determination of scattering events along the detector dimensions 
within its position resolution. In the best cases these events can be observed as a smooth curve 
which reproduces the neutrino disappearance probability. We call this measurement a "neutrino 
oscillometry" . It is worthwhile to note again that the oscillometry is suitable for monoenergetic 
neutrino since, then, one deals with a single oscillation length L32. This is obviously not a case for 
antineutrino, since, in this instance, one extracts only an effective oscillation length. Thus some 
information may be lost due to the folding with the continuous neutrino energy spectrum. 



EXPERIMENTAL APPROACHES TOWARDS NEUTRINO OSCILLOMETRY 

One could propose different experimental methods, already mentioned above, which can be carried 
out for the scan of neutrino scattering events. One suggestion is to install the neutrino source in the 
centre of sphere and to measure the events in the 47r-geometry which provides the highest angular 
efficiency, leadingg to the geometric factor g s = 1 as discussed above. This, however, is only suitable 
for relatively small L32-values because the sphere dimensions are limited, at least very lsarge TPC 
spheres are not presently available. Another method concerns to the use of long detector (cylinder 
type) with installation of the source on the top and/or the bottom of detector. Though angular 
efficiency is much smaller than in the case of sphere, g c <C g s , the oscillation path could be much 
longer. In addition, the neutrino source can be removed each time, if necessary, in order to perform 
the background measurements or to install a new sample for repeated measurements. 
In the calculations of the expected number of events for both these methods one needs to know the 
neutrino intensity of a specific source (see Table [XJ) and the neutrino-electron elastic scattering cross- 
section. The dependence of this cross-section on neutrino energy is not sensitive to the elemental 
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TABLE I: Possible candidates for neutrino box-oscillometry. Only the nuclides with half-lives shorter than 
10 5 y and with the neutrino production intensity higher than 10 12 s^ 1 have been chosen. 



Nuclide 


Ti/2 




E v 


L 23 /2 




weight 


V- 




(keV) 


(keV) 


(m) 


(keV) 


gr 


intensity(s _1 ) 


41 Ca 


10 b y 


421 


417 


208 


260 


400 


10 12 


bi Cr 


28d 


753 


747 


373 


560 


250 


10 18 




2 7v 


232 


226 


110 


106 


4000 


5 x 10" 


71 Ge 


lid 


232 


222 


110 


100 


300 


2 x 10 is 


lU3p d 


17d 


543 


480 


240 


315 


15 


5 x 10 ib 


lU9 Gd 


460d 


214* 


101 


50 


30 


50 


5 x 10 ib 


i39 Ce 


138y 


113* 


74 


37 


20 


1.5 


2 x 10 14 


ib 'Tb 


70y 


60.0(3) 


9.8 


5 


0.4 


5 


2 x 10 14 


lby Ho 


4500y 


« 2.6 


w 0.5;« 0.8 
;« 2.2;^ 2.3 
2.6 


0.2-1.3 


< 0.03 


250 


5 x 10 i2 


iyj pt 


50y 


568.0(3) 


44(70%) 
54(30%) 


22 
27 


6.5 
9 


300 


5 x 10 i4 



composition of the detector's target material and is shown in Fig. [T] for the energy interval up to 
750 keV. 




FIG. 3: The rate ^ (per meter) for Ar at 10 Atm with 1 Kg of 157 Tb (a) and 193 Pt (b) as a function of 
the source-detector distance (in m ). The results shown correspond to sin 2 2^13 = 0.170,0.085 and 0.045 
(decreasing from bottom to top). This rate was obtained for a running period equal to the half life of the 
source. 



Spherical gaseus TPC-detector with the internal monoenergetic neutrino source 

In this spherical chamber the neutrino source is proposed to be situated in the centre of the sphere 
and the electron detector is also placed around the source in the smaller sphere with radius r = 0.5 
m. The sphere volume out of the detector position is filled with a gas (a noble gas such as Ar 
or preferably Xe, which has a higher number of electrons). The recoil electrons are guided by the 
strong electrostatic field towards the Micromegas-detector [lj], The NOSTOS, Time Projection 
Counter, has an advantage in precise position determination (better than 0.1 m) and in detection of 
very low electron recoils (down to a few hundred of eV, that suits to all nuclides from Table U except 
of 163 Ho). 
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We have seen that the number of electron events is given by Eq. [I7l This includes: 

• The elastic neutrino-electron cross section. 

(see Eqs fTS] and [TB")) . The scale is set by the value given in Eq. [3] 

• The number of neutrinos emmitted by 1 Kg of the source. 



1 m s 6.0 

1 Kg A s x 1.66 x 10" 27 ~~ 1 Kg~A s 

where A s is the atomic number of the source and m s is its mass. For any given time t the 

£n2 



^ = T^ , ^ C :^,- 27 = T^-?X^ (28) 



number of neutrinos emitted must be multiplied by the fraction 1 — e~ l ' T with r = -j^t)- 



• The number of electrons present in the target: 
Assuming that we have a gas target under pressure P and temperature To we find: 

P 97 , P Z 300 

n e = Z- — =4.4xl0 27 m -3 29 

kT 10 Atm 18 T y ' 

(Z the atomic number). 

By using this value, as well as the neutrino intensity and the neutrino-electron cross section (see 
Fig- P), we estimated the expected total neutrino scattering events as a function of distance from the 
centre of sphere, for sources 193 Pt and 157 Tb, including the oscillation effect. Here we assumed that 
the sphere radius was equal to r=10 m and the Ar target was under pressure P=10 Atm at room 
temperature. We have considered three values of the parameter sin 2 (2^13), namely sin 2 (2613) — 
0.170, 0.085, 0.045. These values were chosen in accordance with the estimate of sin 2 6i^ — 0.02±0.01 
recently obtained [l?} • 

The obtained results are shown in Figs and [5Jd. In these plots one should pay attention not 
only to the absolute value of the rate but also to the dip due to oscillation. Clearly very accurate 
measuremets are needed. From the length of oscillation point of view the ideal choice is 157 Tb, since 
the full oscillation takes place inside the detector, but 193 Pt could be a good compromise. 



Cylindrical detector with the external monoenergetic neutrino source 

An alternative approach consists in using a long cylinder filled with a liquid target, which can 
provide much longer scanning length. Since the liquid target has higher density than a the gaseous 
one of the spherical TPC, one expects higher event rates. Admittedly, however, for an external target 
installed on the top of the detector, there will be flux losses caused by smaller effective acceptance 
angle, if the whole range of length L is to be explored. This will be espacially true for large L, in 
which case the cone with the vertex at the place of source will be smaller than a half sphere, let 
alone a full sphere. As a matter of fact we have seen that the geometric factor, compared to that of 
a half sphere, is smaller (see Fig. [5]). These losses are in some sense compensated by the possibility 
to scan the long L32-valu.es which are attributed to relatively short-lived nuclides, which results to 
higher neutrino source intensities that can be produced (see Table HJ . 

The use of an external and reproducible source has very big advantage because allows to measure 
the background when source is taken away. This background for the proposed experiment can be 
generated by cosmic and partly by geo-neutrinos. As an example one can estimate the neutrino 
scattering rates for the cylindrical detector LENA [l5| which can be filled by liquid scintillator in 
the tank with the diameter 30 m and height 90 m. The parameters of this detector allow to reach the 
expected position resolution of ca. 60 cm and the energy resolution of 40 keV for 200 keV energy. 
The best cosmic background conditions are reached at 700 keV energy for the LENA-detector. 
Therefore in Table U we include the 51 CY with the neutrino energy 747 keV with a correspondingly 
long oscillation length in comparison to the length of detector (90 m). One may be able to cope 
with this problem by satisfying oneself with only a portion of the oscillation inside the detector. The 
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strong source of 51 Cr, of at least 30 MCi, can be prepared for this case. Figs [5} [8] show estimated 
total rates for the electron- neutrino scattering for sources 51 Cr, 55 Fe, 71 Ge and 103 Pd for the above 
values of sin 2 (28 13 ) . In estimating the value of A the electron density of n e = 2 x 10 29 m -3 for 
LS-LAB will be used. 

In this case the presentation of the results is more complicated, since the total rate, for a given 
neutrino energy, depends on two variables p and C (see Eq. (l2~Tj) ) and the geomtric factor is a bit 

complicated. We will present our results in units of A = N„n e R (see Eq. ([23|) ). Thus only the 

dependence on the neutrino energy and the geometry will be exhibited. 

Naively one might expect to see the oscillation patter in the case of cylindrical geometry directly 
from the raw data in a fashion analogous to that of the sphere. This naive view, however, is not true 
as seen in Fig. [4] The variation of the neutrino flux from point to point outweighs the variation of the 
oscillation amplitude from point to point. Clearly a much more careful analysis should be employed. 
This is indeed achieved by counteracting the effects of the flux variation, i.e. by multipling the event 
rate by (£ 2 + u 2 p 2 ). This way we obtain results as follows: 




C = z/h — » p = r/R — > 

(a) (b) 

FIG. 4: The differential event rate , in units of A (see text), for a cylinder of radius r — 15m and length 
h — 90 m and E v — 226 keV. On the left panel (a) the plots are preseted as a function of £ = z/h for the 
values of p = r/R — (0.05, 0.1, 0.15, 0.2, 0.5, 1.0), to be read from top to bottom. On the right (b) the same 
quatity is presented as a function of p = r/R for the values of £ = z/h = (0.05, 0.1, 0.15, 0.2, 0.5, 1.0). The 
effect of oscillation is not visible. It is included only for orientation purposes, i.e. to demonstrate that a 
more careful analysis is needed. 



• Rates for the source 55 Fe (neutrino energy 226 keV): 

In this case from the data of table fl] and the above value of n e we obtain A = 0.243s -1 = 
7.7 x 10V 1 - 

If we multiply the obtained event rate by the above mentioned factor, (£ 2 + u 2 p 2 ), we get the 
results shown in Fig. [5]- [5] From Fig. [5] we see that the magnitude of the rate depends strogly 
on the radial distance from the axis, but the oscillation period depends only mildly on it. 

• Rates for the source 51 Cr (neutrino energy 747 keV). 
In this case A = 0.577s" 1 = 1.82 x 10 7 y _1 

• Rates for the source 71 Ge (neutrino energy 222 keV). 
In this case A = 1.266s" 1 = 4.0 x 10 7 y _1 

• Rates for the source 103 Pd (neutrino energy 480 keV). 
In this instance A = 0.0282s" 1 = 891y _1 

It is rewarding that, if the data analyzed as above, the oscillation pattern exhibited by a cylindrical 
detector is not very different from that of the spherical detector. It should be kept in mind that the 
exhibited range of C is to guide the eye. The maximum value of £ is, by definition, unity. 
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C = z/h - 

(c) (d) 

FIG. 5: The differential event rate (u 2 p 2 + C 2 )^^j (i n units of A) for a cylinder of Radius R = 15m and 
length h = 90m as a function of £ = for 55 Fe (E„ = 226 keV). The solid, dotted and dashed curves 
correspond to sin 2 2#i3= 0.170, 0.085 and 0.045 respectively. The panels (a), (b), (c), (d) correspond to 
values of p = r/R =0.05, 0.2, 0.6 and 0.8 respectively. 



Analysis of the data of a cylinrical detector in terms of L 

In the previous section we analyzed the data by obtaining the event rate for each point p, £ of 
the detector. One may attempt to present the data in terms of the source detector distance L by 
a suitable transformation from the variables (r, z) to (L, </>) and integrating over the angle tf>. A 
tedious but straightforward algebra yields: 

R-^- = N„n e R^g, v (u, L/R)a(E u , L) = A^g, v (u, L/R)u{E v ,L) (30) 
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FIG. 6: The same as in Fig. [5] for 51 Cr (E v = 747 keV). 



where g av (u, L/R) is a geometric factor that takes care of the variation of the neutrino flux in the 
various positions described by L. It can be cast in the form: 



1. 



g av (u 7 v) = 1 - \/v 2 - l/v, 



< v < 1 
1 < v < 1/u 



(31) 



- Vv 2 - 1)/V, l/u<V< y/l + l/l 



This funcion for u = 1/6 is shown in Fig. On the same figure we show the expected event rate 
for a typical case like 55 Fe. It is clear that to disentangle the oscillation dependence on L one must 
divide the expected rate by the geometric factor. When this is done we obtain the results shown in 
Figs [10] and [IIJ 
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DISCUSSION 



In this article we consider two possible experimental approaches for the neutrino oscillometry: 
the spherical gaseous detector and the liquid scintillator cylindrical one. Both cases have their own 
advantages and disadvantages. 

• The spherical gaseous detector. 
This requires the construction of a sphere with the radius much larger than 10 meters can 
hardly be built. This restricts the possible candidates for the neutrino oscillometry by nuclides 
whose neutrinos energies should not be much higher than 10 keV in accordance with the 
equation [27l The relevant oscillometry curves are shown in Fig. [3] employing an Ar gas. One 
kg of source, installed at the center of the sphere, was employed under the assumption of a time 
measurement longer than the half-life of the source. To obtain the number of events collected 
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FIG. 8: The same as in Fig. [5]for 103 Pd (E v = 480 keV). 



during the reasonable measurement time, the values in Figs [3] must be multiplied by a number 
of disintegrated nuclides during this time. Thus, e.g., for one year measurement with the 0.15 
kg of 157 Tb one can obtain the rate equal approximately to 8 x 10~ 4 events per meter, too 
low value. The similar feature is expected for other relatively weak neutrino emitters listed in 
the bottom part of Table |U Thus we can conclude that the spherical gaseous detector with 
the neutrino source in the center, being very attractive in its idea, cannot provide presently 
the appropriate conditions for neutrino oscillometry. As a matter of fact, the requirement 
of low energy neutrinos brings the long-lived candidates with the low rate of events. The 
small neutrino energy is associated with the tiny neutrino-electron cross-section, and the small 
electron recoil energy requires the use of gaseous environment with the TPC-detection which 
has not a high target electron density. All these disadvantages with an additional demands for 
very long irradiation at the reactor for the production of appropriate nuclides and a very long 
data acquisition time makes this approach very problematic. 
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FIG. 9: In panel (a) we show the geometric factor <; av (l/6, L/R) as a function of L/R for the value u = 
R/h = 1/6. Note the flux factor of 1/2 has explicitly put in in Eq. (|30p and thus it is not included in the 
geometry. In the first region, which is a spherical segment, one sees that the geometric factor is unity like 
the case of the sphere. In panel (b) we show the expected event rate (in units A = 0.243s -1 ) in the case of 
the 55 Fe source as a function of L in meters. It is obvious that both plots follow the same pattern. 




FIG. 10: In panel (a) we show the expected event rate R^r, divided by the geometric factor, in units of 
A = 0.243s -1 expected for the source 55 Fe, while in panel (b) we show the same quantity for 51 Cr in units of 
A = 0.577s -1 . In both panels the solid, dotted and dashed curves correspond to sin 2 29n~ 0.170, 0.085 and 
0.045 respectively. One does not see the full oscillation inside the detector, but the information is adequate 
to extract usefull information on the neutrino oscillation parameters. 



• A long cylindrical detector. 
In this case the neutrino source in a form of a small sphere can be installed at the top/bottom 
of the cylinder at the center of the circular face. Even though the angular efficiency of such 
configuration is much less than for 4-7r acceptance of the spherical detector, there are many 
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advantages which overbalance this lack. If the length of cylinder is rather long (let's presume 
about 100 m) the neutrino emitters with the energies in the hundreds keV region can be used 
for oscillometry. This on one hand increases the neutrino-electron cross-section (see Fig. [T|, 
while on the other drastically increases the neutrino source decay rate because of decreasing 
the source half-life. Then the increased in the energy of the recoiling electrons makes possible 
to use a liquid target, which leads to the increase of the number of target electrons. 

As it can be seen from the top part of Table Q] there are several candidates ( 51 Cr, 55 Fe, 71 Ge, 
and 109 Pd) for whom the neutrino flux is a few order of magnitudes higher than for the 
heavier nuclides. The (v e ,e) cross-section for, e.g. 71 Ge or 55 Fe, are at least two orders of 
magnitude higher than that for 157 Tb, a typical candidate for oscillometric measurements in 
the spherical gaseous TPC. One-two orders of magnitude can be won for electron density in the 
liquid target. Thus, approximately eight orders of magnitude gain for a long cylindrical liquid 
detector is against about of one-two orders of magnitude angular efficiency loss in comparison 
to spherical gaseous case. Meanwhile, there is some disadvantage in the use of cylinder whose 
length is restricted by 100 m or a little bit more. For relevant candidates collected in the U 
with the neutrino energy higher than 220 keV the scan in a full oscillation length (including 
both disappearance and appearance events) is impossible. Meanwhile, as it can be seen from 
Fig. El the oscillation curves can be well recognized (till z/h =1) and the dependence on the 
mixing angle #13 is also visible. 

In order to show the feasibility of neutrino oscillometry for long cylindrical liquid target let's 
comment the curves in Fig. [5] These curves have been calculated for the neutrino energy of 226 
keV in the cylinder with the length h = 90m and with the radius R = 15m filled by the liquid 
scintillator (LS). The scale parameter A = 0.243s -1 appearing in these curves can be obtained with 
assumption that the neutrino source intensity is N v = 5 x 10 17 s -1 (see Table H| and an electron 
density 2 x 10 29 m -3 via Eq. [23l As a factor u 2 p 2 + rj 2 — (r 2 + z 2 )/h 2 it is less than 1 for practically 
full neutrino path in the cylinder, the range of values of the average differential rate, d 2 N/dpd(, 
is from > 10~ 3 s _1 - for the top left curves of Fig. \5\to > 2 x 10 _2 s _1 for the bottom right ones. 
The values of differential accounts strongly depend on the (r, z)-position of event in the cylinder. 
Another evidence for successful oscillometry can be shown from the presentation of the data in the 
d 2 N/dpd( -values. Though the expected position resolution in the cylinder with LS is better [l5| 
than 1 m , we will assume the volume resolution equal to 1 m 3 . Fig. [5] shows an experimental 
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"fingerprint" which is expected from the existence of neutrino oscillations. 

The solar neutrino background depends on the energy region of interest and it becomes small 
above 400 keV.The expected background from the solar events in this region is < 1 day -1 m -3 for an 
LS detector of the size used in our calculations. It is thus useful to obtain some oscillation features 
in this energy region, using the two sources 51 Cr and 103 Pd, which are favorable in this region. The 
oscillation features of these nuclides are shown in Figs [5] and [5] One can obtain at the end of the 
cylinder a difference of 8.7 x 10 _3 i?s _1 for values sin 2 29 13 = 0.170 and 0.045 in the case of 51 Cr. This 
yields 5.8 x 10 m s -1 . Taking into account the radioactive decay of 51 Cr we obtain 1560 events 
per meter for 60 days of measurent. This difference exceeds the 4<r standard statistical uncertainty, 
even for this unfovarable candidate. Note that the estimated rates in the liquid scintillator presented 
in Fig. [5] are higher than expected background mentioned above. 

More impressive results are expected in the case of 55 Fe and 71 Ge, which are characterized by 
the fact that almost the full oscillation takes place inside the detector (see Figs [5] and [7| as well as 
1101 and ITT ]) . We should mention, however, that in this energy one encounters higher solar neutrino 
background [lij. 

Quite high values obtained can be even increased if the data handling is going during several 
months or even weeks. These measurements can be stopped in any time. Source can be moved 
and new source installed. Because of high neutrino intensity the reasonable result can be obtained 
very quickly without the agonizing long-term measurements. This advantage of removable sources 
is indeed very attractive. 

CONCLUSIONS 

The new method for neutrino oscillation measurements is proposed. It is based on the use of 
monoenergetic low energy neutrinos which are released in the atomic electron capture by nuclei. 
High intensity neutrino sources with the neutrino energies less than a few hundreds keV can be 
produced in neutron reactors. They can be installed at the top of the long cylinder filled by liquid 
target material or placed in the center of gaseous sphere placed underground. The measurement 
consists in finding the number of events arising from the electron neutrino scattering on the electrons 
of the target. It was shown that, for some appropriate candidates like ( 55 Fe, 71 Ge, 109 Cd), the so 
called small oscillation length L32 can fall within the dimensions of cylindrical detectors which can 
be constructed in the nearest future. There the oscillations can be determined by disappearance of 
events caused by the change of the neutrino flavor just inside the detector. These events can be 
scanned point-by-point thus providing the oscillometry curve. The calculations for the cylindrical 
detector with the length ft = 90 m and radius R = 15 m showed a complicated net of events 
which can be analyzed by the "oscillation signature" suggested. The most attractive advantage of 
proposed experiment is the possibility to move the neutrino source, to measure the background or to 
implement other physics project and also to replace one source by another. With the neutrino source 
intensity of 30 MCi, which provides high count /background rate, the observation of oscillations can 
be performed during months or even weeks in underground conditions . Since the $13 discovery 
potential for proposed experiments spans the values of sin 2 29i3 < 0.170, the experiments with low 
monoenergetic neutrinos, e.g., with the LENA-detector, could be complementary to investigations 
with the reactor and accelerator neutrinos in Double Chooz, RENO, Daya Bay, T2K and NOvA 
projects planned for the coming decade. The physics interest in the short baseline oscillometry 
measurements is not only the determination of #13. The discovery of the L32 -value from the 
oscillometry experiments with the subsequent comparison with the electron neutrino energy will 
allow to check availability of equation (|27[) , the dynamics of v e — e low energy interactions and, in 
general, the validity of the long baseline oscillation results obtained so far. The present analysis does 
not suffer from the 8-fold degeneracy [l9| that corrupts the measurement of #13 in conventional long- 
baseline experiments. The neutrino energies, needed to check (|27j) . can be accurately determined by 
independent direct mass measurements of nuclides, a capture- partners, in the ion traps fllj ]. 
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